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Abstract 

In this paper, we introduce a family of observables for the dimer model on a bi- 
periodic bipartite planar graph, called pattern density fields. We study the scaling limit 
of these objects for liquid and gaseous Gibbs measures of the dimer model, and prove 
that they converge to a linear combination of a derivative of the Gaussian massless free 
field and an independent white noise. 

1 Introduction 

1.1 The dimer model 

Let G be a graph. A dimer configuration C of a graph G is a subset of edges of G such that 
every vertex of G is incident to exactly one edge of C. The dimer model is a system from 
statistical mechanics, obtained by endowing the set of all possible dimer configurations 
with a probability measure. It has been introduced in the 1930s || to give a model for 
adsorption of diatomic molecules (dimers) on the surface of the crystal, represented by the 
graph G. This model of statistical mechanics is one of the rare models that can be solved 
exactly. For an introduction to the dimer model, see for example Q. 

We will require G to have at least one dimer configuration. Suppose for the moment that 
G is finite. One can define a Boltzmann probability measure on the dimer configurations of 
G as follows : positive weights w e are assigned to the edges e of G and the probability of 
a dimer configuration C is chosen to be proportional to the product of the weights of the 
edges it contains: 

p[e]=^n™e, (!) 

eee 

where the normalizing factor Z = X^erieee^e * s caue d the partition function. 

Kasteleyn showed || that if G is planar, the partition function Z can be expressed 
as the Pfaffian of a weighted adjacency matrix for a well-chosen orientation of the graph. 

"Centrum voor Wiskunde en Informatica - P.O. Box 94079, NL-1090 GB Amsterdam 
cedric .boutillier@cwi .nl 



1 



If moreover the graph is bipartite, Z reduces to the determinant of a certain matrix K, 
called the Kasteleyn operator, a cousin of the adjacency matrix of G, whose lines are 
indexed by white vertices, and columns by black vertices. In particular, it has the property 



that if e = (w,b), |K(w,b)| = w[e}. See [12] to see how to define K. Because of the 



correspondence between determinants of submatrices of a matrix and those of its inverse, 
the correlations are given by determinants of submatrices of K _1 ||. 

When the graph G is infinite, it may have an infinite number of dimer configurations, 
and it is not possible anymore to define directly a Boltzmann measure. For planar bipartite 
Z 2 -periodic graphs, endowed with periodic weights on edges, this notion is replaced by that 
of Gibbs probability measure that has the following properties: 

• it is ergodic under the action of Z 2 by translation, 

• if the dimer configuration is fixed in an annular region, then the random dimer config- 
uration inside and outside the annulus are independent, and the induced probability 
measure inside the annulus is the Boltzmann measure defined above. 

Sheffield proved [|14|] that for given periodic weights, there exists a family of Gibbs 
measures leading to the same Boltzmann measure on finite regions parametrized by the 
two components of an external magnetic field B. We will only consider here the case 
when B = 0, since the other measures of this family can be obtained as measures without 
external field for different weights. 

The fundamental domain G\ of G is supposed to have the same number of white and 
black vertices. The white (resp. black) vertices are labeled, say, from 1 to n: wi, . . . , w n 
(resp. bi, . . . , b n ). If v is a vertex of G, then v x = v + (x, y) denotes the translate of v by 
the lattice vector x = (x, y) S 1? . If (z, w) is in the unit torus T 2 = S 1 x S , a function / 
is said to be (z, it;)-periodic if 

VveG, V(x,y)GZ 2 , f(v x>y ) = z- y w x f(v). 

As in the case of a planar finite region, a infinite Kasteleyn matrix K is defined. Its 
action on functions defined on black vertices returns a function on white vertices 

(K/) w = K ( w > b )/b 

w~b 

If e = (w,b), we will often write K e instead of K(w,b). 

Since K is periodic, its Fourier transform K(z,w), being the restriction of K to the 
(z, u;)-periodic functions, is a n x n matrix. The lines of K(z, w) are indexed by the white 
vertices of the fundamental domain G\ (more precisely, by the (z, w)-periodic function 
taking value 1 at a given white vertex of the fundamental domain, and elsewhere). Its 
columns are indexed by the black vertices of G\. 

The correlations between dimers under the Gibbs measure are expressed in terms of 
determinants whose entries are coefficients of the inverse Kasteleyn operator K -1 given by 
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the usual inverse Fourier Transform formula. Let b and w be respectively a black and a 
white vertex in the fundamental domain, and x = (x, y),x' = (a/, y') be lattice points. The 
coefficient K _1 (b x ,w x ') is given by 

K-(b,.w*) = K-(b_„w) = 

z -(y-y ) w {x-x )Q ifW [z, w) dz dw 
•2 P(z, w) 2iirz 2mw 

where K~ 1 (z,w) is the inverse of the n x n matrix K(z,w), P(z,w) its determinant and 
Qh,vj{ z -,w) is the cofactor associated to K w ^(z, w). 

The probability that edges ei = (w 1 , b 1 ), . . . = (w fc , h k ) appear in the random dimer 
configuration is given by 

F[ ei ,...e k }= (r\K(^,bi)) det [K _1 (b*, w- 7 ')] . (2) 

L=i J 

The way the correlation between edges decays with the distance depends on the number 
of zeros of P(z,w) on the unit torus Three behaviours, or phases, are possible for a 
Gibbs measure. The measure can be: 

• solid, when there are bi-infinite duals path on which edges appear with probability 
or 1. 

• liquid when K _1 (b x , w) decays linearly with |x|, and thus correlations decay polyno- 
mially. A more precise statement is given in lemma B 

• gaseous when K~ 1 (b x ,w) and thus correlations, decay exponentially. 

We will here deal with fluid Gibbs measures, corresponding to the last two cases. 

1.2 Scaling limits of pattern densities 

A dimer configuration on a bipartite planar graph can be interpreted through the height 
function as a discrete surface ||16|| . From this point of view, scaling limits of dimer models 
on planar bipartite graphs have already been the object of several studies: a law of large 
number has been established 0, [ll]] showing that this discrete surface approaches a limit 
shape when the mesh size goes to zero. The fluctuations of the height function around the 
limit shape have also been studied in the case of a graph embedded in a bounded region, 
as well as in the case of an isoradial infinite graph with critical weights Q . The continuous 
limiting object for these fluctuations is the Gaussian Free Field [13, [|]. 



In this paper, we are interested in the scaling limit of dimer models on Z 2 -periodic 
planar graphs but from a different standpoint. Instead of looking at the height function, 
we consider other observables, called pattern density fields. 
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Let G be a Z 2 -periodic planar bipartite graph. A geometric realization of G is an 
application from G to R 2 preserving the Z 2 -periodicity of G: vertices of G are mapped to 
points of M 2 , edges to segments, and Z 2 acts on the image of G by translation. 

Let $ be a geometric realization of G such that the quotient of R 2 by the action of 
Z 2 has area 1. For each scaling factor e > 0, we define the scaled geometric realization 

= e^f and G £ = ^ £ (G) the image of G by the application * £ . 

A pattern IP is a finite set of edges {ei, . . . , e^}, together with a marked vertex. The 
position of the pattern in G £ is given by the coordinates of the image by \E ,£ of this marked 
vertex. If ej goes from white vertex wj to black vertex hj , then the probability of seeing 
such a pattern in a random dimer configuration is given by formula (||) 

P[3>] = fjlKCw^oj JeyK-^w,)]. 

In order to get some information about spatial distribution of patterns, and the way they 
interact with each other, we define for every pattern IP a family of (discrete) random fields 
Ny, called pattern density fields. For a given e, Ny is a random distribution, associating 
to every domain D the number of copies of IP seen in D in a random dimer configuration 
of G £ . More precisely, if u £ is the image by *$> £ of the marked vertex attached to IP, the 
action of Nj, on a smooth test function <p is given by 



m<p) = E <p«)*9* = > ; (3) 



xez 2 xez 2 



where IP X and are the translates by x of IP and u £ . To simplify notations, IP X will also 
represent the indicator function Iy x of the pattern IP X , equal to 1 or whether IP X is in the 
random dimer configuration C or not. 

e 2 E [iVy (<£>)] is a Riemann sum of P[IP](/j, and thus converges to P[!P] f (p(u)\du\ when e 
goes to 0. 

The aim of the paper is to prove the following convergence of the fluctuations of this 
field around its mean value: 

Theorem 1. When e goes to zero, the normalized fluctuation field 

N £ 7 = e{m-nm\) 

converges to a Gaussian field, which is the linear combination of a derivative of the Gaus- 
sian Free Field and an independent white noise. 



See section 1.4 for a more precise statement. Such central limit theorems exist for 



linear statistics of a broad class of determinantal random fields [15] with Hermitian kernels. 
However, since the kernel K _1 appearing in dimer models is not Hermitian, those results 
do not apply. 

Before stating precisely the main results, we recall some basics facts we will need about 
Gaussian fields. 
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1.3 Gaussian Fields 



A Gaussian field is somehow a infinite dimensional generalization of the notion of Gaussian 
vector. See || for an introduction. As in the classical situation with a Gaussian vector, all 
the moments can be expressed in terms of the second moment, by Wick's formula: 

Proposition 1 (Wick's formula). Let W be a Gaussian field. All the moments of W 
are determined by the covariance: Let (p±, . . . , ip n be smooth test functions. Then 



E[W(ipi)---W(v n )) 



if n is odd, 

n/2 

£ l[E[W(ip lk )W(ip jk )] if n is even. 

_ t pairings k=l 



(4) 



In particular, if all the test functions are taken to be equal and E [^((/j) 2 ] = a 2 , then 
we recover the usual formula for the moments of a Gaussian random variable: 

E[»W] = (° „ i !" iSOdd ' 0» 

I [n — 1)!! a if n is even. 

where (n — 1)!! = (n — 1) ■ • ■ 3 ■ 1 is the number of pairings of n elements. 

To prove the convergence in distribution of a sequence of random fields (W n ), in the 
weak sense, to a Gaussian random field W, one has to prove that for all smooth test 
functions (p, W n ((p) converges in distribution to W(tp). Since the distribution of W(<p) is 
determined by its moments, one has just to check that the moments of W n {<p) converge to 
that of W(ip), i.e. that Wick's formula is asymptotically satisfied by (W n ). 



1.4 Statement of the result and outline of the paper 

The main result of the paper is the following Central Limit Theorem when the measure on 
dimer configurations is fluid (i.e liquid or gaseous). 

Theorem 2. Consider the dimer model on a planar bipartite 7? -periodic graph G with a 
generic liquid or gaseous Gibbs measure [x. Let T be a pattern of G and NS, be the random 
field of density fluctuations of pattern y. Then when e goes to 0, Ny converges weakly in 
distribution to a Gaussian random field N;p. 

• If the measure is liquid, then Ny is a linear combination of a directional derivative 
of the massless free field and an independent white noise and its covariance structure 
has the following form 

E pVpOpipV 3>(¥>2)] = ol II d(pi(u)dif2(v)G(u, v )dudv + (3 Lpi(u)tp2(u)du (6) 
where G(u, v) = — log \u — v\ is the Green function on the plane. 
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• If ijl is gaseous, then is a white noise and 



E[Xy(^i)Wg>(^ 2 )] =0 f ^(^^(ujdu (7) 

In other terms, for any choice of (pi, . . . , (p n S C£°(R 2 ), the distribution of the random 
vector (N5,(pi, . . . , Nytp n ) converges to that of the Gaussian vector (Nyy?i, . . . , ]%></?„) whose 
covariance structure is mentioned in the theorem. As the distribution of a Gaussian vector 
is characterized by its moments, it is sufficient to prove the convergence of the moments of 
(Nytpi, . . . , NyLp n ) to those of (Nyipi, . . . , 3\fy<^ n ), given by Wick's formula (||). 

The proof goes in two steps: first we prove the convergence of the second moment of 
the fluctuation field Ny to the covariance of 3%>, and then we prove that Wick's formula is 
satisfied asymptotically. 

As the arguments are different for a liquid and a gaseous measure, the proof of theorem ^ 
is decomposed into three cases. In section ||, we give the proof for a pattern made of a single 
edge in the generic liquid case (when the two zeros of P on the unit torus are distinct) and 
discuss briefly what happens in the non generic case. In section [|, the proof is extended to 
any admissible pattern ( i.e. a pattern appearing with positive probability) in the generic 
liquid case. The situation when the measure is gaseous is discussed in section |j| The 
correlations between different pattern density fields are presented in section |6|, providing a 
generalization of theorem ^. In the last section are given two explicit computations : one 
on the square lattice, and the other on the so-called square-octagon graph. Before entering 
into the details of the proof, we present in section || some properties of the different fluid 
phases of the dimer model. 

2 The fluid phases of the dimer model 

We present here some properties of liquid and gaseous phases of the dimer model on a 



planar bipartite periodic graph G. See [12] for more details. 



2.1 Liquid phase of dimer models 
Map from G* to M 2 . 

When given a liquid measure on dimer configurations of G, there is a natural application 
ty* from the dual graph G* to R 2 , described by the following lemma. This application 
seems to give the good geometry to study liquid dimer models. In particular, when dimer 
weights are critical, ^* coincides with the isoradial embedding of G* described in ][l0|]. 

In the liquid phase, the characteristic polynomial P has two zeros on the unit torus, 
that are complex conjugate of each other, and generically distinct. Let {zq,wq) be one of 
them. 
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Lemma 1. The 1-form 

e = (w,b) i ^ iK v , h (z ,w )Q hv ,(z , w ) (8) 

is a divergence- free flow. Its dual is therefore the gradient of a mapping from G* to M. 2 ~ C . 

This mapping is 1? -periodic and the symmetries of its range are generated by x = 
iz diP(z ,w Q ) and y = iw d 2 P(z , w ). 

Proof. The divergence of the 1-form u : e i— ► iK Vf ] 3 (zo, wo)Qh v/ (zQ, wo) at some black vertex 
b is given by 

(divw)(b) = ^ w 'b(^o, wo)Qbw'(^o,wo) 

w'~b 

= «(Q(z ,^o) • ^(zo,wo)) b)b = iP(2o,wo) = 0. 

Similarly, one can check that the divergence of this flow is at every white vertex w. Thus, 
since G is planar, there exists an application : G* — > C such that uj = d^. 

The fact that ^f* is periodic is a consequence of the fact that G and the Kasteleyn 
operator are both periodic. There exist two complex numbers x and y such that for every 
/ G G* and every (x,y) E Z 2 , the difference between the image of f x>y and that of / itself 
is given by 

-**(/) = ** + Jtf. (9) 

Define a and (3 to be the partial derivatives of P at the root (zq,Wq) with respect to the 
first and the second variable respectively 

a = diP(z , w ), (3 = d 2 P(z , w ). (10) 

Let us prove now that the numbers x and y are given by iz^a and iwof3 respectively. Let j x 
and 7^ be the path delimiting respectively the lower horizontal boundary and the leftmost 
vertical boundary of the fundamental domain G\ (figure ||). The complex number x equals 
the sum of the complex numbers ±u;(e) over all edges e crossing ^y x (see figure 0). These 
are exactly the edges of the fundamental domain whose weights have been multiplied by 
z ±l in the Fourier transform of the Kasteleyn operator K(z, w), and the sign of the power 
of z is the same as that in front of u>(e). Noticing that 

z if m = 1 

-z~ x ifm = -l (11) 
if m = 

we can write 

x= ^2 *( 2; o9iKwb(2o,'Wo))Qbw(^o,'Wo) 

w,bGGi 

= iz tr (diK(zo, w ) ■ Q(z ,w )) ■ 



z ^- z m = mz m 
oz 
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Figure 1: A piece of a planar bipartite periodic graph G . The shaded zone represents its fundamental 
domain Gi, delimited by the two paths y x and ~f v . 



On the other hand, since the characteristic polynomial P(z, w) is the determinant of 
K(z, w), then 

n 

iz a = iz diP(z ,w ) = iz ^ det(K hl (z , w ), . . . ,d 1 K hj (z ,w ), . • . , K h Jz ,w )) 

i=i 

(12) 

where K^. (z, w) is the j'th column of the matrix K(z,w). Expanding each determinant 
with respect to the column containing derivatives, we get 

n n 

iz a = izo^2^2 di^w^. (z ,w ) x 
j=i k=i 

Cof Wfcbj {K hl (z , w ), . . . , dxK^ZQ, w ), . . . , K hn (z , w j) ■ (13) 

Since the cofactor Cof Wfc b 3 {K bi , . . . , diK^^ (zo,wq), . . . , K\> n (zq,wq)) does not depend on 
the jth column, we can replace it by K\ 3j (zo, wo). This cofactor by definition is Qbj\v k (zo,wq) 
Thus, 

n n 

iz a = iz diK Wkbj (z , wo)Q hjWk (z , w ) 

j=l k=l 

= iz tr(d 1 K(z , w Q ) ■ Q(z , w )) = x. (14) 
The same argument applied to 7 y gives the formula 

i WQ P = y. (15) 
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□ □ 

In what follows, to construct the application from G* to M. 2 , we will choose between 
the root (zo,wq) on the unit torus for which the frame (x, y) is direct. This is equivalent 
to requiring that 

Im (^)>0. (16) 

To get a geometric realization of G from Vl/*, just pick a point in each dual face of 
^*(G*) in a periodic way. Note that sometimes, as in the case of isoradial embeddings [1C] 



the obtained realization is not a plane, but can be thought as a globally flat manifold with 
conic singularities, obtained by gluing together the images of faces of G* . 

Asymptotics of K _1 (b,w). 

The coefficients of K" 1 decay linearly. More precisely, if b and w are in the same fundamen- 
tal domain, and b x>y is a translate of b by (x,y), then we have the following asymptotics 
for K^b^, w): 



Lemma 2. Let (zo,wo) the root of P on the unit torus satisfying (jig). Then the asymptotic 
expression for the coefficients of K" 1 are given by 

U-KU X r> ( ZQ~ y WQ X Q h ^(zQ,W )\ ( 1 

K L (b x ^,w) = -Re — - -. '— r— } + 



m{azox + (3woy) j \\x\ 2 + \y\ 2 
WoQb,w(^o,w ; 

7r(xx + yy) / \\x\ A + \y 



Re| V^^^ , +Q( _^ , (17) 



The proof of this lemma is given in |l^ 

Note that the geometry of the map from G* to M? appears in this analytical result. 
The denominator up to a factor it is the vector separating the image of the fundamental 
domains of b XtV and w under this mapping. Moreover, if b and w are the ends of an edge 
e then (|l7|) can be rewritten as 

K(w,b) K-\b x , y ,w)=Re( Z f W !_ iG * s )+o( , |2 ] ■ |2 ) (18) 

\ir[xx. + yy) J \\x\ 2 + \y\ 2 J 

and thus correlations between edges decay polynomially with the distance. 
2.2 Gaseous phase 

In this phase, the characteristic polynomial P has no zeros on the unit torus. If b and w 
are a black and a white vertex of G\, the fraction Qbw/P is analytic on the unit torus and 
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its Fourier coefficients K 1 (b x ,w) decay exponentially with x: 

3d,C 2 >0, VxeZ 2 , Vb,weGi, |K- l (bx,w)| < Cie" C2|x| . 
Hence the correlations also decay exponentially. 

2.3 Density fields and partition function 

In this paragraph, the calculations are purely formal and try to give some heuristics on the 
information one could get from these pattern density fields. 

Consider a dimer model for which we assign to a configuration C the weight w(G). Let 
Zq be the partition function of this model 

z = J>(e) (19) 

e 

We now perturb the partition function by modifying locally the configuration weights. Let 
<p be a smooth test function. Fix a pattern 7 and a e > 0. We multiply every weight 
w(C) by a factor e teip<yU ^ whenever there is a copy of a given pattern 5> located at x. Up to 
a multiplicative constant exp(teP[y] X^^^x))) the new P ar titi° n function for the model 
with these new weights is 

e e 

This can be generalized to a perturbation involving several patterns. Formally, the suc- 
cessive derivatives of ZtjZ§ at t = are the moments of the random variable N^{(p) with 
respect to the unperturbed probability measure. 



d " Z; =E 1^;)' 
t=o 



dt k Z { 



Thus, an information on the moments of gives an infinitesimal information on the 
perturbed partition function that could be hopefully integrated to construct probability 
measures corresponding to the modified weights. 



3 The liquid case: edge densities 

We first concentrate on the proof of a simple particular case of theorem ^. We suppose 
that the probability measure on dimer configurations of G is a generic liquid measure, and 
consider the fluctuations of the density random field iV| for a pattern consisting in a single 
edge e = (w,b). The precise statement we prove in this section is the following: 
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Theorem 3. The random field N% converges weakly in distribution, as e goes to 0, to a 
Gaussian random field N e with covariance 

E[X e ^i)We(^ 2 )] = - jj d e *cp 1 (u)G(u,v)d e *<p 2 (v)dudv + A J \p x (u) V2 {u)du (20) 

for a certain A > 0, and where e* is the vector representing the dual edge of e in the 
geometric realization ty* of G* presented in lemma |^. 

Note that the particular geometry introduced by the application \t is particularly well- 
adapted to the problem. 

As we discussed in the previous section, we first prove the convergence of the second 
moment, and then that of higher moments. In this section, as we are interested in copies 
of an edge e = (w,b), the only vertices we will deal with are most of the time translates 
of w and b. To simplify notations, we will write K _1 (x — x') instead of K -1 (b x , w x /) and 
K e will stand for K(w,b). 

3.1 Convergence of the second moment 

The second moment ((fi,(p 2 ) l— ► IE N^ipx)N^((p 2 ) of iV| is a continuous bilinear positive 

form on C£°(IR 2 ). We prove that this bilinear form converges to a non degenerate bilinear 
form, that will define the covariance structure for the limit Gaussian field N e . 

Proposition 2. There exists a non-negative constant A such that 



limE 

e^O 



N^(ipi)N^(ip 2 ) =-ll d e *ipi(ui)G{ui,U2)d e *ip2(u2)\du 1 \\du 2 \ 



IT 



A / ¥>i(u)¥> 2 (u2)|du| (21) 



where G(u, v) = — Jj- log \u — v\ is the Green function on the plane. 

Before going into the proof of this proposition, we give some interpretation of the 
expression for the covariance. The right hand side can be physically interpreted as the 
energy of interaction between two magnetic dipoles with moment density ^ie* and </?2 e * > 
plus a term of interaction at very short range. 

Suppose that there is an excess of edges e in the random dimer configuration, in some 
region D\. These dimers behave collectively as a magnetic dipole: their presence influences 
the rest of the dimer configuration as if a magnetic field created by a dipole with a density 
e* over D\ was applied to the system: edges whose dual has an orientation closed to that 
of the magnetic field at that point is more likely to appear. 
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Proof. Using the invariance by translation of the Kasteleyn operator and hence of the 
correlations, we rewrite the second moment as a convolution of two distributions, applied 
to a test function. 



= £ 2 E ^i(M x )^2(u x ')E[(e x -e)(e x , - e)] (22) 

= e 2 Y ^i(ux)<^2(wx') E [( e - e)(e x - x ' - e)] (23) 
x,x'ez 2 

= (24) 
where e = P [e] and the two distributions ip\ and F £ are defined by 

ip\ =e 2 ^2i Pl (u x )<5 u e F e = E [(e - e) (e x - e)] 5 U ^ . 

xy x 

The distribution (pf converges weakly to (pi when e goes to zero. We will now prove the 
convergence of F £ to some distribution F, what will ensure that tp\ * F £ converges weakly 
to tp\ * F 6 , since the support of tp\ is contained in the fixed compact supp(^i), and hence 
that {(p\ * F £ , ip2) converges. 

Let tp be a smooth test function with compact support. Let us prove the convergence 

of 

(F £ , V) = E E t( e " ^)( e x " e)] Y>«) = E Cov ( e > e *M*4)- (25) 

X X 

At first sight, (F £ ,ifj) looks vaguely like a Riemann sum of a particular function. The 
problem is that, due to the asymptotics of K _1 , the function would behave as 1/u 2 , which is 
not integrable in the vicinity of 0. Therefore, we decompose the sum on x in the definition 
of F £ depending on whether the norm of x is larger than M = |l/ e J ° r n °t, that is if u x 
is in B = {— i(azos + 0wot) ; (s, t) £ [—1, l] 2 } or not. 

<F £ ,V>= E V'(n x )Cov(e,e x )+ Y ^«)Cov(e,e x ) 

|x|>M |x|<M 

= Y, V>«)Cov(e,e x )+ £ (V'«)-V'(0))Cov(e,ex)+V(0) E Cov(e,e x ). (26) 

|x|>M |x|<M |x|<M 

The fact we subtracted and added ^(0) in the second sum removed the non integrable 
singularity at 0. The following two lemmas state the convergence of the three sums. 



E 
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Lemma 3. 



lim V V«)Cov(e,e x ) = -L f 4>{z)Re(—^ 



|x|>|l/e| 



1 



lim (^«) " ^(°)) Cov ( e ' e ^ = ^2 ~ ^(°)) Re - l du l 

£ " |x[<Ll/eJ 



(27) 
(28) 



Moreover the sum of the two previous limits can be rewritten as 



- d e *i/j(u)d e *G(u,0)\du\ +Ai^(0). 

/or some constant A\. 

Proof. If x ^ (0,0), the covariance between edges e = e( ,o) and e x = is given by 

P[e and e x ] - P[e] 2 = K^det 
-Kl K-^K-^-x) 



Cov(e, e, 



K- J (0) K- X (x) 
K- J (-x) K- x (0) 



(KeK-^O)) 5 



Using asymptotics of K 1 for large x, we get the following asymptotic expression for 
the covariance between two distinct edges 



Cov(e, e> 



r Re 



z y w$ie* 



7T 



-lii 



(I- 



2vr 2 



Re — + 



22/ ""o 2:r ) 



+ o - 



Since the second term is oscillating, it will not contribute to the limit. The sum in the 
left hand side of (p7D, modulo the oscillating terms, can be interpreted as the integral of a 
piecewise constant function approximating 

-_L^(u)Re 

As the approximating functions are bounded uniformly in e by an integrable function, 
and converge almost everywhere, then by Lebesgue theorem, the first sum converges to 



(29) 
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In the same way, the sum in the left hand side of ( p8[ ) is the integral of a piecewise 
constant function approximating 

-^i(V'(n)-V(0))Re H(u) 

and for the same reasons, it converges to 

We rewrite the sum of the limit using Green's formula inside and outside of !B, noticing 
that <^:Re (^-) is the second derivative of the Green function G(u, 0) = — ^ log \u\ along 
the vector e*. 

-d 2 e ,G(u,0)iP(z)\du\ = 



- (f i>(u)d e *G(u,0)(n in ,e*)da - - [ d e * 

7T JdT, K JR 2 \3 



4>(u)d e *G(u,0)\du\, (31) 



-d 2 e ,G(u,0)(Hu)-m)\M = 
n if 

-<f (^(«)-^(0))a e *G(u,0)<fr tmt ,e*)d(r-i I 8 e *iJj(u)d e *G( U ,0)\du\. (32) 

where fti n and n OU £ are the unit normal vector fields on d'B pointing respectively inwards 
and outwards. The two integrals on 23 and M? \ H combine to give an integral over M?. 
The two integrals on d'B cancel out partially. It remains only 



-V>(0) i> d e *G(u,0){n ext ,e*)da = A^O). 
J (923 

This establishes the convergence of the two first sums , which completes the proof of the 
lemma. □ □ 

We now prove the convergence of the third sum. 

Lemma 4. ^| X |<M Cov(e, e x ) converges when M goes to infinity to a limit A^. 

Proof. The sum of the covariances is given in terms of K and K -1 by 

Cov(e,e x )=P[e](l-P[e])+ ^ Cov(e,e x ) 

\x\<M 0<]x|<M 

= K e K-\0)-Kl Y, .K-^xJK-^-x) (33) 

|x|<M 
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K *(x) is by definition the xth Fourier coefficient of the function / = Q/P defined on 
the unit torus T 2 . As P has simple zeros, / is in L 1 (T 2 ). The convolution 

du dv 



(f*f)(z,w) 



f(u,v)f(zu,wv) 



(34) 



2mu 2ittv 

is also in L X (T 2 ) and its xth Fourier coefficient is exactly K _1 (x)K _1 (— x). Establishing the 
convergence of the sum is now a problem of pointwise convergence of a Fourier series. If 
/ had been continuous at (z,w) = (1, 1), then the Fourier series would have converged to 
/(l, 1). The problem is that the function f(u,v) 2 is not integrable and thus, the function 
/ * / is not defined when z and w are both equal to 1. However, / * / is smooth in a 
punctured neighborhood of (1,1), has directional limits when (z,w) converges to (1,1), 
varying continuously with the direction. We can then prove an analogue in two dimensions 
of Dirichlet's theorem for / * / to show that the Fourier series at (z,w) = (1, 1) when M 
goes to infinity, converges to a mixture of the directional limits of / * /. More precisely, if 
t = &rg(w)/ arg(z) and £(t) the limit of / * / when (z,w) goes to (1, 1) with t fixed, then 



M- 



1 + t 



1 - t 



dt 

T 



(35) 



i r+oo 

lim y K-\x)K- 1 (-x) = ^ £(t)log 

l^+oo ^ IT 2 J „o 

|x|<M x 

And thus 5Z|x|<M Cov(e — e)(e„ 

We now come back to the proof of the convergence of the distribution F £ . The three 
sums in (p6|) defining (F £ ,^) converge and the sum of the limits is 



e) converges to a limit that we denote by A<i- □ □ 



(F,V> = - / d e *^(u)d e *G(u,Q)\du\ +A^(0) 

TT 



(36) 



where A = A\ + Ai- Thus, when e goes to 0, F £ converges to the distribution F defined 
by the formula above, and hence <pf * F £ to (pi * F. Denoting by <9 (u) and d M respectively 
the operator of partial differentiation with respect to the variable u and v, and noticing 
that, since G(u,v) = G(u — v), we have: 

d {u) G{u,v) = -d (v) G(u,v), 

we get the following expression for the limit covariance structure 



limE 



(Jv^iXiv^) 
i 



((pi * F, ip 2 ) 
d e * (pi (u)d^G(u, v)p 2 (v)\du\ \dv\ + A 

d e *(pi(u)d^G(u,v)(p 2 (v)\du\\dv\ + A 

R 2 xR 2 

d e *(pi(u)G(u,v)d e *(p2(v)\du\\dv\ + A 



(pi(u)(p 2 (u)\du\ 

(pi(u)(p 2 (u)\du\ 
(pi{u)(p 2 (u)\du\ 



(37) 



Thus the covariance of N £ converges to the expression given in proposition 



□ □ 



15 



3.2 Convergence of higher moments 



We now prove the convergence of the moments of order > 3 of iVJ to those of the Gaussian 
field ]NT e . 

Proposition 3. For every n > 3, The nth moment of iV| converges to that of N e when e 
goes to zero. In other words, for every <fi, - ■ ■ ,<f n S C£°(M 2 ), 



limE 

£->0 



Since N e is Gaussian, it is sufficient to show that in the limit, the moments of satisfy 



JV^i ■ ■ ■ Nfy n 



is a symmetric ra-linear 



Wick's formula. Moreover, as (<pi, • • • , ip n ) i— > E 

form, we just have to prove proposition [| when all the tpi are equal to some test function 
tp, the general case being obtained by polarization. The previous proposition reduces then 
to showing that 



Proposition 4. 



limE 







[ (n - 1)!! E [(^ e V) 
In this section, we are thus interested in the limit of 



2W2 



if n is odd 
if n is even 



E 



^1 )*** 1^71 

A first step in the proof is to study the convergence of a related quantity E^(ip), defined 
by a sum of the same general term as for E 
to distinct points: 



e" 2^ ^i u . 

XI— Xn 
distinct 



xi - 



(N^xfj) n , but with a set of indices (xj) restricted 
^«)E[(e X] -e)...(e x „-e)]. (39) 



3.2.1 Convergence of S^(^) 

To prove the convergence of E^(ip), we have to understand the asymptotic behavior of the 
correlations between distinct edges, when they are far from each other. A simple expression 
is given by Kenyon in |J to compute these correlations using a unique determinant. 

Lemma 5 ([^]). Let e\ = (wj, bi), • • • , e Ti 

is given by 

K-^b^Wj) 

E[(ei -ei)---(e n - e n )] = I J|K(w j ,b j 



w n ,b n ) be distinct edges. Their correlation 



U =1 



det 

l<i,j<n 



K-^bj.Wi) 
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This formula allows us to give an explicit expression for S^(V') m terms of the operators 
K and K -1 . Since the matrix in lemma |5| has zeros on the diagonal, only permutations with 
no fixed point will contribute to the determinant expressed as a sum over the symmetric 
group. Let & n be the set of such permutations. Every permutation a G 6 n is decomposed 
as a product of disjoint cycles ji' — 'yp. The supports of these cycles form a partition 
(Ti)f =1 of {1, ...,n}, whose parts Ti have cardinal at least 2. The terms coming from 
permutations leading to the same partition are put together and we get: 



xi---x n 

distinct 



J]^«.)(e x .-e) 



3=1 



ik ™ e n^K) det 



xi ■ ■ -x n j=\ 

distinct 



K — 1 (b,-,»j) 



Xl-'-X„ ,i jf = l 

distinct <Jfc0n 

/ 



En e 

{r}f =1 i=l 7 cycle . 

supp(7)=r ; \ 



I- 1 



sgn( 7 )^l E n^K fc ) K e K "(x 70ft) -x ife ) +o(l 



distinct 



\ 



(40) 



The error term o(l) comes from the fact that in the second line, we allow two Xj whose 
indices are in different components of (T[) to be equal. 

We now examine the convergence of a term in brackets, associated to a cycle 7 



sgn( 7 )^l £ J] ^.JKeK- 1 ^) - x jk ) 
x ji '"' ' x j| 7 | fc=i 

distinct 

According to subsection [O], we know that when 7 is a transposition the corresponding 
term converges. When the length of 7 is at least 3, we have the following lemma. 

Lemma 6. For any cycle 7 of length m > 3, and any ip £ Q^(R 2 ), we have 

m 



distinct 



1 



Re 



dui\ ■ ■ ■ \du m \ (41) 



2 m l^m J (R2)m ^ (u 7(1) - ui) • • • (« 7 ( m ) - U m ) , 
Proof. From the behavior of K -1 at long range, we get an asymptotic expression for the 
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product 

m m 

nK e K- 1 (x 7(j) -x j ) = nRe 

3=1 3=1 



:: + gmall 



o 



7r(it 3 



•7(3) 



terms 



(42) 



-Rc 



(ie*) 1 



2 m-l 7r m I („| u e ) . . . ( u e u e ) 

\ v x 7(l) x l ' v x 7(m) X m /^ 



oscillating 
terms 



(43) 



The oscillating part of this asymptotic expansion once summed over xi , . . . , x m will not 
contribute to the limit. The sum of the leading term multiplied by e m ^(u^ i ) • • • ip(u^ m ) 
can be interpreted as the integral of a piecewise constant function, converging almost 
everywhere to 



f-^tp(ui) ■ ■ ■ ip(u m )Re Yl 



.J|[ '"7O') U 3 

As all the functions are dominated by a constant times the integrable function 

\ip(ui) ■ ■■tp(u m )\ 



(44) 



|(« 7 (i) - - (n 7(m) -u m )\ ' 
the convergence follows from Lebesgue theorem. 



□ 



□ 



Once proven the convergence of all these terms, we can combine their limit to get the 
limit of E^(ip). When summing over all cycles with a given support T = {ji, . . . , j m } of 
cardinality m > 3, we get the following limit : 



lim > 

£^0 ^ 



( 



7 cycle 



e m £ f[^ jk )K e K-\ Xj{jk) -xj +o(l) 



V 



X il '"' > x jm fc=l 



supp(7)=r \ distinct 
1 



2"!— lyj-m 



(R 2 ) r 



E 11 



, 7 cycle JET 

\supp7=r 



U l{j) U 3 



i| " 



/ 



which equals zero according to the following lemma : 
Lemma 7. Let m > 3, and u\, . . . ,u n be distinct complex numbers. Then 



m 1 

E U^- 1 



= 



7G6 m i=l 
m— cycles 



|du r 



(45) 
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Proof. Denote by / the function of u%, . . . , u m defined by the left hand side of (45). When 
m is odd, the result is obvious, since 7 and 7 -1 give opposite contributions. For a general m, 
since m-cycles form a conjugation class in the group of permutations 6 m , the function / is a 
rational fraction invariant under permutation of the variables u±, . . . , u m . The denominator 
of this fraction is the Vandermonde V = Yli^iui—Uj), and the numerator is of lower degree 
than V. Since V is antisymmetric under permutation, the denominator has to be as well. 
But the only antisymmetric polynomial of lower degree than V is 0. □ □ 



In the limit of equation fl40| ) will contribute only partitions whose all components have 
cardinality 2. If n is odd, {l,...,n} cannot be partionned into parts of two elements, 
and the limit H n (?/>) of is zero. When n is even, there are (n — 1)!! such partitions, 

corresponding each to a pairing. The limit E n (ip) of E^(V0 is then 

n/2 / 

S «W= E II LV E ^« lfc )^« Jfe )Cov(e Xifc ,e Xj) 

pairings fc=l \ x ifc ^x jfc 

{h,jl},— ,{i n /2iJn/2} 

= (n-l)H (S 2 (^)r /2 

For the moment, we discussed the limit of a restricted sum on distinct edges, which is 
not exactly the expression for the nth moment. We will now deal with the case when some 
edges can coincide. 

3.2.2 Proof of proposition ^] 

In the expression of the nth moment 

e n Y, ^i)--^«)E[(e xl -e)---(e Xn -e)] 

Xl,- ,Xn 

the correlations are not given by lemma |^ as soon as at least two edges coincide. To 
understand the behavior of this expression, we must be able to express in terms of K and 
K _1 correlations of the form 



E 



(e Xl -e) fc i...(e Xp -e)M (46) 



when ei, . . . , e p are distinct, and k\, . . . , k p > 1. Using the fact that the indicator function 
of an edge e satisfies e J = e for j > 1, Newton formula yields 




= a e k (e-e) + [3l 
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where and /3f are deterministic, depending only on e and k. Since all the edges we 
consider are translates one from another, and thus have the same probability, we will 
simply denote these coefficients by ct k and (3 k . Note some particular values of at- and (5k 
that will be useful later 

at = l, 0x = O, (5 2 = e(l-e). 
Correlation (^) can be rewritten with these notations as 



E 



(e Xl - e) kl ■ • • (e» 



e) kp 



E [{a kl (ei - e) + (3 kl ) • • • (a kp (e p - e) + /3 fep )] 



E ^ E 

Jc{i,-, P } 



n(e,-e) 



(47) 



where 5 j = a fcj . and/3 j = • 

In equation (|38|), E (iV|'0) n J is expressed as a sum over all edges. We want now 

to rewrite this sum as a sum over distinct edges, using partitions of {1, . . . , re}. Such a 
partition {A;}^ =1 is associated naturally to every re-tuple of lattice points (xi,...,x n ): 
each component A; of this partition is an equivalence class for the relation 



l ~ J Xj 



Denoting by re; the cardinal of A;, we rewrite equation (|38| ) as 



E 



E E « 

{A,}f . xi,-,xp 
1 Ji=1 distinct 



[J^'«.)(e x ,-e) 



i=i 



e " E E E E 



{A,}f =1 JC{l,...rf»} 



distinct 



e 'n^ 



(48) 



where £' means the sum over (xj) distinct, but also distinct from values of any Xj,j 6 J. 
Denote by g the number of Aj reduced to a single element. As (5\ = 0, /?j is zero unless 
J contains the indices of these Aj. Thus, the cardinal of a subset J giving a non-zero 
contribution must be at least q. For such a J, the last sum over (x/Wj in ( f48|) is a 
Riemann sum, and therefore is 0(s~ 2 ( p ~\ ,J ^). 

Furthermore, the sum over (xj)j g j can be expressed by polarization in terms of Hfyi, 

and is therefore a 0(e~l J l). Since |J| > q and 



n 



E n i - 1 + ~ i) = 2 p - ^ 



i=i 
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the contribution of J to ( f48| ) is at most 0{e n 2 P+I J 1) which will be negligible in the limit 
except when | J\ = q and n\ = 2 for all I ^ J. For such J and (A/,), we have 



aj = 1, Pj = (e(l - §)) 



(n-g)/2 



(49) 



Thus the only partitions that will contribute to the limit are "partial pairings" , matching 
by pairs 2m = (n — q) elements of {1, . . . , n}. For a fixed m, there are 

n\ 



n 

2m 



(2m - 1)!! 



2 m m!(n - 2m)! 

such partitions, all giving the same contribution. Summing over m we get 



E 



[n/2\ 

E 

771=0 



n 



(2m - Et, 



, )niJ <- _„_ 2m (V)(e 2 e(l-e)^V 2 «)) + 0(e). 

The Riemann sum e 2 e(l — e) V 2 ( u x) converges to 

e(l - e) J iP 2 (u)du = E [H e (^) 2 ] - 3 2 (# 
If n is odd, so is n — 2m. In this case, limH^_ 2m (?/>) = 0, and therefore 



limE 



0. 



If n is even, 



(V) = H n _ 2m (^) = (n - 2m - 1)!! ~ 2 (V>) n/2 - m 



(n — 2m)! 



2™/ 2 - m (n/2 - m)! 



S 2 W 



n/2—m 



(50) 

(51) 
(52) 



Therefore, the limit of E 



limE 

£^0 



n/2 

E 

m=0 
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is given by 
(n-2m)! E 2 (^) n/2 ~ m 



2 m m\(n-2m)\ 2 n l 2 -' m (nj2 - m)\ 

n/2 



e(l-e) / V 2 W|du| 



= (n - 1)!! £ ( n/2 ) 3 2 (^/ 2 -™ (e(l - e) / ^ 2 (u)|du| 

m=0 V 771 / V •/ 

/ /■ \ n/2 

= (n - 1)!! f 3 2 ty) + e(l - e) J ^ 2 (u)\du\) 

= (n-l)!!(E [^ e (^) 2 ]) n/2 

what is exactly what we wanted to prove. This therefore ends the proof of theorem ^ for 
a pattern made of one edge and a generic liquid Gibbs measure. 
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3.3 A word about the non generic case 

When the two roots of the characteristic polynomial P(z, w) on the unit torus coincide, 
the measure is still liquid, and the correlations between edges at distance r still decay like 
r~ 2 . However, since zq and wq are real, the leading term is the asymptotics of K" 1 is not 
oscillating anymore, what will induce a "resonance phenomenon" in the system. 

The two first sums in equation (|2^) defining the distribution F 6 appearing in the study 
of the convergence of the second moment still have a finite limit when e goes to zero. On 
the contrary, due to the resonance, the third sum 

^2 Cov(e,e x ) 

M<[l/eJ 

in this case diverges. More precisely, this sum is 0(log(l/e)). Therefore the second moment 
diverges. However, one can prove that (log(l/e)) _1 iV| converges weakly in distribution to 
a white noise. We will not show it here. 

4 The liquid case: general pattern densities 

This section is devoted to the proof of an analogue of theorem [3| for multi-edged patterns. 
Let Ny the density fluctuation field of a pattern 7 for a generic liquid Gibbs measure. 

Theorem 4. When e goes to zero, Ny converges weakly in distribution to the Gaussian 
field whose covariant structure is given by 

E[3<y((pi)'Ny((p2)] = - if dy*(pi(ui)G(ui,U2)d v *(p 2 {u2)\dui\\du2\ 

+ A ipi(ui)if 2 (u)\du\ 

where the vector !P* and the nonnegative constant A depend only on the Gibbs measure and 
the pattern 7. 

The scheme of the proof is very similar to that of theorem j| for edge density fluctuations. 
The problem is that there is no simple analogue of lemma y for correlations between non 
overlapping patterns, that is why the proof needs a little more combinatorial work in this 
case. 

After having introduced the different notations required to deal easily with these pat- 
terns, we prove the theorem, following the structure of the proof given in the last section, 
and explaining in details only parts that are specific to patterns made of more than one 
edge. 
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4.1 Notations 



Let 7 be a pattern containing k distinct edges e 1 = (w 1 , b 1 ), . . . , e k = (w k ,h k ). The 
probability of such a pattern to appear in a random dimer configuration is 



3> = P pp] = K e i ■ ■ ■ K efe det 



K-^bSw 1 ) 
K-^b*^ 1 ) 



K-^b 1 ^*) 
K _1 (b fc , w fc ) 



(53) 



More generally, the probability to see n non-overlapping copies CPi , . . . , T n of T obtained 
respectively by translation of a lattice vector xi , . . . , x n is given up to a constant by a 
determinant of matrix nk x nk defined by blocks 



P y n ] = (K e i---K efc ) n det 



' All 




Al n 








. A n l 




4 



(54) 



where the entries of the block Ajj are coefficients of K 1 between black vertices of (Pj and 
white vertices of Tj 



Au= : •.. : . (55) 

The matrix A// does not depend on /. We denote by E this matrix, whose determi- 
nant is used to compute 7 = P[!P]. We suppose that the pattern appears with positive 
probability T > 0, and in particular, E is invertible. Defining Bjj as the product E~ 1 Au 
and B as the whole block matrix (Bu), we can rewrite the joint probability of (Pi, . . . , 7 n 
as 





h 


^1,2 




B\ : n 


(?) n det 


-E>2,1 


















-E>n,l 




Bn,n—1 


h 



(56) 



Instead of using a single integer i to denote the line (resp. the column) of an entry 
in such a matrix defined by blocks, it will be more convenient to use a couple of integers 
(I, a), where / is the index of the block line (resp. of the block column) and a is the 
relative position in the Ith block line (resp. block column). The relation between the two 
sets of indices is simply 

i = I(k - 1) + a. 
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If the coordinates Xj are all distinct but some patterns partially overlap, define 
\ Ui=i ^ for all i € {1, . . . , n}. We have then 



?1 ■ • • ?n 



Up to a relabeling of the patterns, we can assume that none of the CPj is empty. Thus the 
joint probability of these patterns is also given by the determinant of a matrix defined by 
blocks of size | V\ \ + ■ ■ ■ + | T n \ . 

4.2 Asymptotics of correlations 

The following lemma gives asymptotic correlations between distant patterns. 

Lemma 8. Let (xj) = xi, . . . ,x n be n distinct lattice points. The correlations between the 
patterns y xi , . . . , J > Xn can be rewritten as 



E [(3> Xl -?)■•■ (0> Xn -?)]= J2 II H v{(*i)) I + 

see n \7 c v cle °f s 
where the functions H~ f have the following asymptotic behaviour 



small 
terms 



2M-i7rl7 



Re 



tr((E- x g)H 



+ 



with Q = (Q havr f3(z ,w )) p =v and satisfy 

\H 7 (( Xj ))\<s^ 



it 

esupp7 u, x 7(j) 



C 



oscillating 
terms 



rijesupp7 \ u i{j) u j\ 

for every u\, . . . , u n in a e -neighborhood of u Xl , . . . , ii Xn . The error term o(l) is uniformly 
bounded in xi,...,x n and goes to zero when the distance between the patterns goes to 
infinity. 

Proof. We first derive the asymptotic expression for the correlations when the patterns 
are far from each other. When I large enough for every i ^ j, the patterns are 

disjoint and expression ( |56[ ) for correlations can be used. Expanding the products in the 
expectation, we get 



E [(? X1 -?)••■ (? Xn -9)]= (-5T~ p l£ 



31,-,3-p 



Cc{l,...,n} 





[ In 


£l2-Bl2 




&lriBln 


det 


^21^21 


In 






















In 



(57) 
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where the non diagonal bloc (I, J) is either Bjj or depending on whether (I, J) belongs 
to C x C or not. 

Expressing each determinant as a sum over the symmetric group & n k and gathering 
the terms coming from the same permutation, one can notice that the contributions of 
permutations fixing a whole block are vanishing, due to the alternating sign in the sum over 
C. The permutations contributing to the correlations are those whose support intersects 
each block. Therefore, we have 

E [(3> X1 -?)•■■ (0> x „ - ?)] = P [0>] n Yl s ^ II %.«).*((/,«)) 

fixing no block 

The main contribution to this sum is given by the "special" permutations, the support 
of which intersects each block exactly once. Let a be such a permutation. The n non 
fixed elements are (1, ot\), . . . (n, a n ). For every / £ {l,...,n}, the non fixed element 
(I, oei) is sent to a((I,ai)) = (S(I), 05(7), where S £ S n is a fixed-point free permutation, 
having the same signature as a. There are k n special permutations a leading to the same 
S, corresponding to the different possible choices of the non fixed points in each block 

Qt\ , . . . , OL n . 

The contribution of the other permutations will be negligible because of the extra e 
coming from additional entries of K -1 between non-fixed points. 



E 



3=1 



9) 



[9] n ]T S gn(a)l[B M(Tait 



a)) 



fixing no block 



I. a 

k 



■pr £ sgn(s) nE( B w)) 



a I, a S(I) 



see n 



\I=1 a I= l 



+ ! ma11 (58) 
terms 



The main term in equation ( p8[ ) has an expression in terms of traces of products of 
block matrices Bjj 

n k 

n e ( B mn) aias(I) = n tr ( B h,i 2 ■ ■ ■ b Wi ) (59) 

•f=l«/ =1 7=(Ji,...,J p ) 

cycle of S 

Let us have a look to a particular trace tr ( Bj 1: i 2 ■ ■ ■ Bj j^. Recall that Bjj is the 
product of E" 1 whose entries will be denoted by e a p and Ajj whose entries (Au) a p are 
the coefficient of K -1 taken between black vertex 'a' of pattern / and white vertex '/?' of 
pattern J. The asymptotics of {Au) a p when patterns are far away from each other are 
given by lemma 0: 



X/' "XjJ 



yj+yi „,,+ x j- x i 



-Re 







Qap(zo,W ) 



7T 



U 



x,/ 



U 



x/ 



+ ! ma11 . (60) 
terms 
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To simplify notations, we introduce 



„ -yj+yi „,.+xj-xi 

UJ, 



cu = p "r — (6i) 



and write Q a p instead of Qy } c, v/ /3(zo,wo). The trace tr (Bj 1: i 2 • • • Bj p j^ can therefore be 
rewritten as 



tr (B h;h ■ ■ ■ B IpJl ) = ^ ( E lA hh)a ia 2---( E 1 Al p l 1 ) ap a 1 
ai,...,a k 

„ , ^ \ small 

terms 



«t V / 



ft,- A 



As in the case of edge densities, there will be only two non oscillating terms in the 
expansion of this product of real parts: those for which the phases contained in the Cu 
compensate exactly. 

With the convention that p + 1 = 1, one has 



tv(B hh ---B IpIl ) = 

/_ £ \p v 

= Yl ( 2W II + QPjaj+iCljlj+i) + 

ai,...,a k \ / 7= i 



small 
terms 



ft,- A 



E(— e) p / t~t _ rr \ oscillating 

11 f ' • 11 ' - /,<^ <;,.,. •>/./,.• + termg 
ai,-,a k \7=1 7=1 

The product of the (i k i k+1 is equal to 



and we can rewrite the trace of -B/ 1) / 2 • • • Bj ^ as 

, /o o \ (-£) p Pr tv{{E-^QY) oscillating 

tr (B, 1>A • • • B IpJl ) - + terms (63) 

giving the asymptotics for if 7 . 

When the patterns are not disjoint anymore, then a similar analysis can be done, in 
defining new patterns as the connected components of "?i U ■ ■ ■ U 3V The bound can be 
extended to this case. □ □ 
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4.3 Convergence of the second moment 
Proposition 5. 



limE 

e^O 



Ny(ipi)Ny{ip 2 ) = -// d v *(pi(ui)G(ui,u 2 )d v *cp2(u 2 )\dui\\du2\ 



7T 



+ A (fi(u)cp 2 (u)\du\ (64) 

JR 2 



The proof of the convergence of the second moment goes exactly as that of section ||. 
The second moment of Nj, can be expressed as a convolution of two distributions applied 
to a test function 



E 



(A^i)(iV^ 2 ) = (<p £ *F s ,ip 2 ) (65) 



with the same definitions as before for y>\ and F e . 



e 



e 2 E mok- - o f £ = Y, Cov ^ ^w- - o ( 66 ) 



(fi converges weakly to (pi . The convergence of F e to a distribution F is proven exactly 
in the same way as in section ^. The only difficulty that could appear is the analogue of 
lemma |] proving the convergence of ^ x Cov(T, T x ). Cov(T, T x ) is a linear combination of 
products of diverse values of K~\ If we interpret these products as Fourier coefficients of 
a convolution of functions, then the convergence becomes more obvious: we saw in section 
H that the function whose Fourier coefficients are the product of two K -1 is not defined 
at (1,1), but has directional limits when (z,w) converges to (1,1), and the sum of the 
Fourier coefficient was an average of these directional limits. When more than two K _1 
are involved, the function is even continuous thanks to the multiple convolutions, and the 
Fourier series converges at (z,w) = 1. 

F £ converges then to the distribution F acting on a test function ip as 

{F,il}) = ~[[ il>(u)Be tT ^ ® - du + A0(O) (67) 



7T 



The complex number representing the vector !P* along which are taken the derivatives 
is a square root of tr(E~ 1 QE~ 1 Q), and after application of Green's formula, we get the 
expression of proposition |5| for (ipi * F, ip 2 ). 
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4.4 Convergence of higher moments 
Proposition 6. Let n > 3, and ipi,...,(p n £ e£°(R 2 ) 



limE 



i%>i)---JV$(¥>„) 







j/n is odd 



n/2 



pairings 1=1 



n is even 



(68) 



Proof. As usual, it is sufficient to study the case where all the ipi are equal to some fixed 
smooth test function tp. The nth moment is the given by 



E 



(my 



E E 



.1=1 



(69) 



We know from lemma the asymptotics of the general term of the sum 



2£ 2|7l 



E II ^7^Re ( tr((^Q)H) J] 



5e<3„ 7 cycle of S 



oscillating small 
+ 

terms terms 



E IK 2 ) 

CO? '=i 



-e 



2\ |r,| 



2tt 



Rc 



v 



7 cycle jgr ; x 7(j) x : 

supp7=r ; 



small oscillating 
+ + 

terms terms 



where the (Ti) are partitions of {!,..., n} whose components have size at least 2. The 



expression we obtained is very close to that of equation (42). From this point, the same 
arguments as for edges yield the proof of the proposition. □ □ 



5 The gaseous case 

In a gaseous phase, K" 1 decays exponentially. There exist two constants C\ and C2 such 
that 
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Vx,x'e z 2 



K-^Voo) 



The precise statement of theorem ^ in this particular context for a pattern consisting 
in a single edge e is the following: 

Theorem 5. The random field JVJ converges weakly in distribution to a white noise of 
a 2 if 



amplitude J ^logK , where 



3"= J J log P(z,w) 



dz dw 



(70) 



2iirz 2iirw 

is the free energy per fundamental domain of the dimer model. 

The proof is exposed in the first two subsections, and the case of a more complex 



pattern is briefly discussed in subsection 5.3 



5.1 Convergence of the second moment 
Proposition 7. 



oo/ n 2\ 



limE 

e^O 



d 2 ? 



(pi(z)<p 2 (z)\dz\. 



Proof. As in section ||, the covariance E N £ ((pi)N £ S(ip 2 
ipf * F £ applied to the test function (p 2 , where 



d log K| 

is a convolution of distributions 



ifl = s 2 («x)*(- -<), F £ = Y1 Cov ( e > e x)^(- " < 

xy x 

(pf converges weakly to ip%. We have now to prove that the distribution F e converges 
/ard 
will be 



toward the distribution F = 9 ^^2 If it is the case, then the limit of the second moment 



(ipi *F,ip 2 ) 



d log K| 



d 2 ? 



d log K £ 



ipi(u)<p 2 (u)\du\. 



If x ^ (0, 0), the covariance between edges e and e x is given by 

Cov(e,e x ) =E[(e 0i0 -e)(e x -e)] = -K2K- 1 (x)K- 1 (-x 



(71) 



(72) 



Let ip be a smooth test function with compact support, and TV a large integer. We 
decompose the sum over x in the expression of {F £ , ip) depending on whether the norm of 
x is larger than TV or not. 
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<F £ ,V>= Y, V>«)Cov(e,e x )+ ^ (V>(t£) - 0(0)) Cov(e, e x ) + V(0) £ Cov(e,e x ). 

|x|>JV |x|<iV |x|<JV 

In the first sum, there are at most 0(e~ 2 ) terms since the support of is bounded. As 
|x| > N, each term in this sum is bounded by some constant times e~ 2C ' 2N . Therefore the 
whole sum is a 0(e~ 2 e~ 2C2N ). In the second sum, since |x| < N, the distance between u%. 
and is less than eN. As tp is smooth , (0(w x ) — ip(0)) is 0(sN). Since there are 0(N 2 ) 
terms in the second sum, it is a 0(N 3 e). If we choose for instance N of order e -1 / 4 , these 
two sums converge to zero, when e goes to zero. 

The third sum is absolutely convergent, the limit of (F £ ,if)} is 

0(0) Y Cov(e,e x ), 

meaning that F £ converges weakly to X^xez^ Cov(e,e x )(5 and the limit of the second mo- 
ment is proportional to the L 2 scalar product. The coefficient of proportionality 

Y Cov(e,e x ) 

can be rewritten in terms of the polynomials P(z, w) and Qbw(z, w) = Q e (z, w). 



Y Cov(e, e x ) = E [(e — e)(e — e)] + ^ E i^Ofi " e) (e x - e)] 
xez 2 x^(o,o) 

= P[e](l-P[e])-K 2 Y K-^xJK-^-x) 

x^(0,0) 

= K e K- 1 (0)-K 2 ^K- 1 (x)K- 1 (-x) 

X 

K e Q e (z,w) KlQ e (z,w) 2 dz dw 



j2 P(z,w) P(z,w) 2 2mz2mw 



Since Q e (z,w) = -rM-P(z,w), we have finally 



<9K e 

dz dw d 2 3 



Y, Cov(e,e„) = K^K^ J f log(F( 2 , M )) 

xGZ 2 



2iriz 2-kzw d log K 2 
□ □ 
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5.2 Higher moments 

Proposition 8. The Wick formula is verified in the limit. 

( 



limE 

e^O 



(n-l)HE (iV, 



n/2 



if n is odd, 
if n is even. 



Proof. As in the case of a generic liquid measure, the proof begins with the study of the 
restricted nth moment £ n (ip) defined by equation (|39[). Lemma || yields an asymptotic 
expression with the same structure as in the generic liquid case (|40|): 



{r}f =1 1=1 7 cycle V 

supp(7)=r( 

l7l v 

E U^<J K ' 1 ^)- x ^+°( 1 n ( 73 ) 

x ji >"' ,x i| 7 | k=l 
distinct 

The contributions of cycles of length greater than 3 vanish in the limit, as an application 
of the following lemma: 

Lemma 9. 

m 

Vm>3, lime™ £ H^.)K-\x j+1 - Xj ) = 

xi,— ,x m j=l 

Proof. Define x' = xi and for j > 2 = Xj — Xj_i . The sum can be rewritten using 
these new notations and invariance by translation of operator K -1 as 



e m E^«') E n^«'+^+... + xPK- 1 (x^) K-^-x'-x^.-.-x^) 

As ip is a continuous function on a compact set, it is bounded. The sum on x' has 
0(e~ 2 ) terms. K- 1 (b_ x ._ x /.„_ x / ,w) is bounded independently from the x^-'s. As K 1 
decays exponentially, the sum on x' 2 , . . . ,x' n is bounded. The whole sum is thus a 0(e m ~ 2 ), 
which goes to zero when e goes to zero as soon as m > 3. □ □ 

Thus E^(V0 converges to E n (^>) = (n — 1)!! ^(ip)- The end of the proof deals with 
collisions between edges, which is identical to what has been done for proposition |], leading 
to the result. □ □ 
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5.3 Patterns in gaseous phase 

Combining the notations and the techniques introduced in section Q to deal with correla- 
tions between patterns, and following the steps of the proof of theorem [|, one can prove 
the following 

Theorem 6. Let IP be a pattern in a dimer model endowed with a gaseous Gibbs measure. 
The random field NS, of density fluctuation of pattern 7 converges weakly in distribution 
to a white noise. 

The proof is omitted here. 



6 Correlations between density fields 

In the previous sections, only fluctuations of the density field associated to one fixed pat- 
tern was considered. One can ask what happens for correlations between density fields 
associated to different patterns. To what extent a high density of some pattern in a given 
region of the plane has an influence on the density of another pattern in an other region ? 

This question is answered by the following theorem generalizing the results of the 
previous sections. 

Theorem 7. Consider a dimer model with a generic liquid Gibbs measure /j,. 

• Let 3>x and T2 be two patterns. The bilinear form E Ny (-)N^, (•) on C£°(IR 2 ) con- 
verges when e goes to zero to a bilinear form E [N^ (-)N;p 2 (•)] . 

• If (j, is a generic liquid Gibbs measure, there exists a constant Ay^ 2 such that for 
every test functions (p\ and ip2, 

1 [Nfr (<£i)N:p 2 O2)] = ~ 1 1 d^<pi(u)G(u,v)dy*<p 2 (v)\du\\dv\ 



+ A 



/ (fii(u)(p 2 (u)\du\ (74) 
Jm. 2 



If /j, is gaseous, there exists a constant A<y x <$ 2 such that for every test functions ip\ 
and tp2, 

E [K Vl ((p 1 )'N'p 2 ((p2)] = A Vl y 2 / (fi(u)(p 2 (u)\du\ 



Let . . . ,y„ be patterns (not necessarily distinct). When e goes to zero, the mul- 



tilinear form E 



iVy (•)••• Ny (•) converges. The limit E fNy^-) • • • N;p n (-)] is given 
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by Wick formula: for every test functions <pi, . . . , <p n , 





n/2 

n 

k pairings k=l 



if n is odd 
if n is even 



The coefficients Ay^ 2 are in general not known in a closed form. However some re- 
lations between them can be found. For example, let v be a vertex of the graph, and 
denote by ei, . . . , e m the edges incident with v. The complex numbers e|, . . . , sum to 
zero since they represent the edges of the dual face v*. Even more, for every e > 0, the 
linear combination N% + • • • + iV| m is identically zero. Indeed, for every x € Z 2 , there is 
exactly one edge incident with v x in the random dimer configuration. Therefore the sum 
of indicator functions (ei) x + ■ ■ ■ + (e m ) x is always equal to 1. This relation between the 
random fields iV|. at a microscopic level yields relations between the different coefficients 



A eiBj . Precisely, 



EE 

i=i j=i 



A 



e » e ) 



(75) 



7 Examples 

The theorems in the previous sections state a convergence of density fluctuation in the 
scaling limit to a linear combination of a derivative of the massless free field and a white 
noise. However, they do not give an explicit form for the white noise amplitude. In this 
section we present some cases for which a closed expression for the white noise amplitude 
can be provided in terms of the weights on edges. The first case is the dimer model on 
the graph Z 2 with periodic weights a, b, c, d around white vertices. The second case is the 
dimer model on the square-octagon graph. 



7.1 Dimer densities on Z 2 

The graph we consider here in the graph Z 2 with a bipartite coloring of its vertices. Weights 
are assigned to edges according to their directions: a, b, c, d counterclockwise around white 
vertices and clockwise around black vertices. If none of the weights is greater than the sum 
of the others, the corresponding dimer model is critical |T(|: the graph can be embedded 
in the plane such that all the faces of the graph as well as those of the dual graph are 
inscribed in circles of a given radius (see figure ^) . The Gibbs measure with no magnetic 
field on dimer configurations is liquid. The dual faces are similar to the cyclic quadrilateral 
with sides a,b,c and d. The area of such a quadrilateral is 

Area = - y/(-a + b + c + d)(a - b + c + d)(a + b - c + d)(a + b + c - d) 
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Figure 2: On the left, a piece of the geometric realization (isoradial embedding in this case) of Z 2 fixed 
by the weights a, b, c, d assigned to edges. The region delimited by the tick dotted contour is a fundamental 
domain, and the coloring-preserving symmetries are generated by x and y. On the right, a white face and 
the quantities (angles and side lengths) related to it. 



and the radius R of its circumscribed circle is denned by the relation 

(ab + cd)(ac + bd)(ad + cb) 



(-a + b + c + d){a - b + c + d){a + b - c + d){a + b + c - d)' 



The fact we chose the fundamental domain to have area 1 leads to the following ex- 
pression for the complex numbers representing the dual edges in the embedding of 1?: 



in 



b* 



ib 



-%CZq 



idzr, 



V2 Area 



v2 Area Area Area 

where (zq,Wq) is the root of the the characteristic polynomial 

j-» z \ b cz 

R{z, w) = a H V dz 

w w 

on the unit torus, with the additional constraint that Im(zo) > 0. 
The following theorem is a particular case of theorem |?]. 

Theorem 8. Let GFF a Gaussian free field in the plane and W an independent white 
noise with unit variance. The vector-valued random field 



N £ c 
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converges weakly in distribution to the vector- valued Gaussian Field 



MA 




fd a A 




1 










w 




WJ 



GFF 



abed 



8ttR 2 Area 



-1 
+1 

v-v 



w. 



Proof. The only point we have to explain is the computation of the white noise amplitude. 
This coefficient can be identified from the limit of the second moment of N% . The proof of 
the convergence of the second moment, as before, goes through the proof of the convergence 
of the distribution F e defined by (|25| ) . The sum over x in the definition of F e is decomposed 
into two parts depending on whether k x belongs to some neighborhood of or not. In 
the general case, we considered the neighborhood B = {i(azos + /3zot) ; (s,t) € [—1, l] 2 }. 
However, here we will take an infinite strip 

8 = {i(az s + 0w o t) ; (s,t) € K x [-1, 1]} . 

The condition on x = (x,y) corresponding to u x £ S is x G Z and \y\ < M where 
M = [1/eJ. We have 



(F s ,^) = V>«)Cov(e,e x )+ £ (^«) - V(0)) Cov(e,e x ) + ^(0) ]T Cov(e,e x ). 

|y|>M |j/|<Af 

The inverse Kasteleyn operator in this case is given by 



\y\<M 



K- 1 (x,y) = K- 1 (b x ,w ) 



z y w 



dz dw 



j2 a + b/w — cz/w + dz 2mz 2iirw 



The fact we have an infinite strip allows us to make use of one dimensional Fourier 
series in the x-direction to compute the third sum. Indeed, K^ 1 (x,y) is the xth Fourier 
coefficient of the function f y (w) defined by 



fy(w) 



y-y 



dz 



i a + b/w — cz/w + dz 2mz 



Hence, for a fixed y, 



^K- 1 (x,y)K- 1 (-x,-y)= I f y (w)f. y (w) 



dw 
2mw 



For y ^ 0, the functions f y and f_ y have disjoint support, therefore the sum above is zero. 
When y = 0, the sum is equal to 



y>-.(s,0)K-.(-*,0)=^(w)^=/'. 



dw 



s±bw\ > i 2m{a + b/w) 2 w 

cw — d I 
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The third sum equals 
V>(0) Cov(e,e x )=V(0) ( 0^(0) -a 2 £ ^ K _1 (2;,y)K _1 (— a:, —y) 

\y\<Mxe1 



x&L 
\y\<M 



adz 



I a+hw |>i (a + b/w)2inw 



a 2 dz 



sd±m\ > i (a + bjw) 2 2mw 



^(0) 



— a6Im(w;o) 
7r|a + 6/7^0 1 2 



The other term coming from the application of Green formula can also be computed 
and turns out to be equal to 



-V(o) 



1 



-oMm(i/)o) \ 



2Area tt 



a + 



w ' 



Hence the variance of the white noise appearing in the limit of N% is given, after some 
calculations, by 

1 — a&Im(u>o) a&Im(iOo)^ abed 



tt \a + b/wo\ 2 
A similar computation for E 



2Area / 87ri? 2 Area 

leads to the same expression, with a 



negative sign. As the expression of the coefficient is invariant under cyclic permutation of 
a, 6, c, d, the amplitudes for the other pairs are easily deduced. □ □ 

When d = 0, the dimer model on Z 2 is equivalent to that on the honeycomb lattice with 
periodic weights a, b, c. One can notice that in this case, the amplitude of the white noise 
vanishes. The interaction between dimers on the honeycomb lattice is purely electrostatic. 
We conjecture that it is true only for that particular model. 



7.2 Dimer densities on the square-octagon graph 

The square-octagon graph is a Z 2 -periodic graph whose fundamental domain is presented 
in figure |3[ It contains four white and four black vertices. When every edge is assigned a 
weight equal to 1, the characteristic polynomial is given by 



P(z,w) 



det 



1 







1 

1110 

z 1 w 
-10 11 



1 1 

w 

z w 



When the magnetic field is zero, or weak enough, the dimer model is in a gaseous 
phase. The fluctuations of the density field of an edge can be therefore computed by 
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Figure 3: A portion of the square-octagon graph. Edges are oriented from white end to black end, except 
those whose orientation is represented on the figure. 



taking derivatives of the free energy of the system with respect to the weights, as explained 
in section ||. 

To compute for instance the amplitude of the white noise in the limit density of the 
edges (wi , bi ) , we assign to these edges a weight e a and to the others a weight equal to 
1, and compute the second derivative of the free energy 3~ associated to this model with 
respect to a. For these new weights, the characteristic polynomial is now 



P a (z,w) = det 



e a -i" 

W 2 

1110 
z 1 w 
-10 11 



4 + e a - e a z - - - e a w 



1 

z 



1 

w 



If a is small enough, P a (z,w) has no zeros on the unit torus, and the dimer model is in 
still a gaseous phase in absence of magnetic field. In every point of this gaseous phase, the 
free energy is constant and given by 



J= JJ \ogP a (z,w) 



dz dw 
2m z 2inw 



log(4 + e a ) + 



log I 1 - 

[0,2^ V 4 + e a 



e a +i + e _ i0 + e a +i4> + e _ i0 



d£d0 
2vr 2vr 



Performing the change of variables a = , f3 = and moving the contour of integration 
over a from [0, 2tt] to [—ia/2, —ia/2 + 2ir] using analyticity and periodicity in a, we finally 
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get an expression of J in terms of an absolutely convergent series: 

4e a/2 



5F a = log(4 + e a ) + 



log 1 



[0,2ti 



4 + e a 



cos (a) cos(/3) 



da d/3 
2vr 2vr 



= log(4 + e«) - £ t 



1 / 4e a / 2 



k=l 



k \ 4 + e a 



oo / 

= log(4 + e ")-£ 

fc=l \ 



e°/ 2 
4 + e° 



2k 



(*0V 



(76) 
(77) 
(78) 



3~ a can be expressed as the value of a certain generalized hypergeometric function. The 
Taylor expansion up to order 2, involving the complete elliptic integrals K and E 



a + 



2vr 



gives information on the statistics of the copies of edge (wi,h>i). The constant coefficient 
is the free energy of the initial model, the coefficient of a is the probability of (wi,bi) 

and the coefficient of a 2 /2 gives the amplitude of the white noise describing the scaling 
limit of the fluctuations of the number of edges (wi,b>i) 



limE 



(^wx.bxjV) 5 



2vr 



u) 2 dii. 



Similarly, one can compute the probability of seeing an edge of a square, for example 
(w2,bi), and the amplitude of the white noise: 

P[(w 2 , bl )] = i + ^iM), 

2A-(1) 



limE 

£^0 



57T 



u) 2 du. 



The fact we see elliptic functions showing up is related to the fact that the spectral 
curve {P a (z,w) = 0} in this case is genus-1 algebraic curve. 
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